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HACK’S LAW IN A DRAINAGE NETWORK MODEL: 

A BROWNIAN WEB APPROACH 

By Rahul Roy, Kumarjit Saha and Anish Sarkar 

Indian Statistical Institute 

Hack [Studies of longitudinal stream profiles in Virginia and Mary¬ 
land (1957). Report], while studying the drainage system in the Shenan¬ 
doah valley and the adjacent mountains of Virginia, observed a power 
law relation l ~ a 0 ' 6 between the length l of a stream from its source 
to a divide and the area a of the basin that collects the precipita¬ 
tion contributing to the stream as tributaries. We study the tributary 
structure of Howard’s drainage network model of headward growth 
and branching studied by Gangopadhyay, Roy and Sarkar [Ann. Appl. 

Probab. 14 (2004) 1242-1266]. We show that the exponent of Hack’s 
law is 2/3 for Howard’s model. Our study is based on a scaling of the 
process whereby the limit of the watershed area of a stream is area 
of a Brownian excursion process. To obtain this, we define a dual of 
the model and show that under diffusive scaling, both the original 
network and its dual converge jointly to the standard Brownian web 
and its dual. 

1. Introduction. River basin geomorphology is a very old subject of 
study initiated by Horton [14]. Hack [13], studying the drainage system 
in the Shenandoah valley and the adjacent mountains of Virginia, observed 
a power law relation 

(1) l ~ a 0,6 

between the length l of a stream from its source to a divide and the area 
of the basin a that collects the precipitation contributing to the stream as 
tributaries. Hack also corroborated this power law through the data gathered 
by Langbein [17] of nearly 400 different streams in the northeastern United 
States. This empirical relation (1) is widely accepted nowadays albeit with 
a different exponent (see Gray [12], Muller [22]) and is called Hack’s law. 
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Mandelbrot [21] mentions Hack’s law to strengthen his contention that “if 
all rivers as well as their basins are mutually similar, the fractal length-area 
argument predicts {river's length) l / D is proportional to ( basin's area ) 1 / 2 ” 
where D > 1 is the fractal dimension of the river. In this connection, it is 
worth remarking that the Hurst exponent in fractional Brownian motion 
and in time series analysis arose from the study of the Nile basin by Hurst 
[15] where he proposed the relation l± = I 9 ' 9 as that governing the width, Zj_, 
and the length, F, of the smallest rectangular region containing the drainage 
system. 

Various statistical models of drainage networks have been proposed (see 
Rodriguez-Iturbe and Rinaldo [26] for a detailed survey). In this paper, we 
study the tributary structure of a two-dimensional drainage network called 
the Howard’s model of headward growth and branching (see Rodriguez- 
Iturbe and Rinaldo [26]). Our study is based on a scaling of the process 
and we obtain the watershed area of a stream as the area of a Brownian 
excursion process. This gives a statistical explanation of Hack’s law and 
justifies the remark of Giacometti et al. [11]: “ From the results, we suggest 
that a statistical framework referring to the scaling invariance of the entire 
basin structure should be used in the interpretation of Hack's law." 

We first present an informal description of the model: suppose that the 
vertices of the d-dimensional lattice Z rf are open or closed with probability p 
(0 < p < 1) and 1 — p, respectively, independently of all other vertices. Each 
open vertex u E Z d represents a water source and connects to a unique open 
vertex v E Z d . These edges represent the channels through which water can 
flow. The connecting vertex v is chosen so that the dth coordinate of v is 
one more than that of u and v has the minimum L\ distance from u. In 
case of nonuniqueness of such a vertex, we choose one of the closest open 
vertices with equal probability, independently of everything else. 

Let V denote the set of open vertices and h{ u) denote the uniquely chosen 
vertex to which u connects, as described above. Set (u, h{ u)) as the edge 
(channel) connecting u and h{ u). From the construction, it follows that the 
random graph, Q = (V,E) with edge set E := {(u, h(u)} : u E V}, does not 
contain any circuit. This model has been studied by Gangopadhyay, Roy 
and Sarkar [10] and the following results were obtained. 

Theorem 1.1. Let0<p<l. 

(i) For d = 2 and d = 3, Q consists of one single tree almost surely, and 
for d > 4, Q is a forest consisting of infinitely many disjoint trees almost 
surely. 

(ii) For any d>2, the graph Q contains no bi-infinite path almost surely. 

In this paper, we consider only d = 2. Before proceeding further, we 
present a formal description for d = 2 which will be used later. Fix 0 < p < 1 
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and let {-B u : u = (u(l),u(2)) € Z 2 } be an i.i.d. collection of Bernoulli ran¬ 
dom variables with success probability p. Set V = {u £ Z 2 : B u = 1}. Let 
{[/ u : u e Z 2 } be another i.i.d. collection of random variables, independent 
of the collection of random variables {B u : u € Z 2 }, taking values in the set 
{1,-1}, with P (U u = 1) = P (J7 U = —1) = 1/2. For a vertex (x,t) € Z 2 , we 
consider ko = min{| k\ : k € Z,B^ x+kt+1 \ = 1}. Clearly, ko is almost surely 
finite. Now, we define 

( (x + k 0 ,t + 1) € V, if (x — ko,t + 1) V, 

h(x,t) := < (x - k 0 ,t + 1) € V, if (x + k 0 ,t + 1) ^ V, 

^ (x + U( x ^ko,t + 1) G V, if (x ± ko,t + 1) G V. 

For any k > 0, let 

h k+1 (x,t) := h(h k (x,t)) with h°(x,t ) := (x,t), 

C k ( x ,t ) := ( k )^ V: h k (y,t — k) = (x,t)}, if (x,t) € V, 

\ 0, otherwise, 

C(x,t) := |J C k (x,t). 

k> 0 

Here, h k (x,t ) represents the “fcth generation progeny” of (x,t), the sets 
C k (x,t ) and C(x,t ) denote, respectively, the set of kth generation ancestors 
and the set of all ancestors of (x, t); C(x, t) = 0 if (x, t ) ^ V. In the terminol¬ 
ogy of drainage network, C(x,t) represents the region of precipitation, the 
water from which is channelled through the open point (x,t) (see Figure 1). 
From Theorem l.l(ii), we have that C(x,t ) is finite almost surely. 

Now, we define 


L(x, t ) := infjA: > 0 : C k (x, t ) = 0}, 

as the “length of the channel,” which as earlier is finite almost surely. We 
observe that for any (x, t) € Z x Z, L(x, t) > 0 and the distribution of L(x, t ) 
does not depend upon (x,t). Our first result is about the length of the 
channel. We remark here that Newman, Ravishankar and Sun [23] has a 
similar result in a set-up which allows crossing of paths. 

Theorem 1.2. We have 

lim \fnF(L(0, 0) > n) = — ^—=, 
n—too Toy 71 " 

where ^ := 7 g(p) = ■ 
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OM) 



Fig. 1. The bold vertices on the line y = t — 3 constitute the set Cs(x,t) and all the bold 
vertices together constitute the cluster C(x,t). 


Next, we define 

rk ( x ^ , = | max{;u :(u,t-k)e C k (x,t)}, 
lk(x t ) , = | min {u : (u,t — k) £ C k (x,t)}, 
D k (x,t) := r k (x,t ) - l k (x,t). 


if 0 < k < L(x,t), (x,t) € V, 
otherwise, 

if 0 < k < L(x, t), (x, t ) € V, 
otherwise, 


The quantity D k (x,t ) denotes the width of the set of all fcth generation 
ancestors of (x,t). We define the width process (s) and the cluster 

process Kn' t \s ) for s > 0 as follows: for k = 0,1,... and k/n< s < (k+l)/n, 


( 2 ) 


7oV n 7oV^ 

7ov n 7o V n 


where 70 > 0 is as in the statement of Theorem 1.2. In other words, Dn’^ (s) 
[resp., is defined £> fc (x, t)/{p/Qy/n) [resp., #C k (x,t)/( 70 y^)] at 

time points s = k/n and, at other time points defined by linear interpo¬ 
lation. The distributions of both Dn' J> and K are independent of (x,t). 

To describe our results, we need to introduce two processes, Brownian 
meander and Brownian excursion, studied by Durrett, Iglehart and Miller 
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[7]. Let {VF(s) : s > 0} be a standard Brownian motion with W( 0) = 0. Let 
ri := sup{s < 1 : W(s) = 0} and 72 := inf{s > 1 : W{s) = 0}. Note that ti < 1 
and T 2 > 1 almost surely. The standard Brownian meander, W + (s), and the 
standard Brownian excursion, VFq~(s), are given by 


(3) 

w + (s) 

\W(ti + s(1 - Tl))| 

VI - d 

s G [0,1], 

(4) 

W+(s) 

_ \w (ti + s(r 2 - Ti))\ 

V r 2 - n 

s G [0,1] 


Both of these processes are a continuous nonhomogeneous Markov process 
(see Durrett and Iglehart [6] and references therein). Further, W + ( 0) = 0 
and, for x > 0, P(VF + (1) < x) = 1 — exp(—x 2 /2), that is, W + (l ) follows a 
Rayleigh distribution. 

We also need some random variables obtained as functionals of these two 
processes. In particular, let 

Iq := / IT+ (t) dt and Mq~ := max{fFo~(i) : t € [0,1]}. 

Jo 

Louchard and Janson [20] showed that, as x —>• oo, the distribution function 
and the density are, respectively, given by 

/? r 7 r ) /CK 

P(/ ( ] h > x) ~ —^xexp (—6x 2 ) and fj+(x)~ x 2 exp (— 6x 2 ). 

The random variable Mq is continuous, having a strictly positive density 
on (0, oo) (see Durrett and Iglehart [6]) and for x > 0, 

OO 

P(Mg" < x) = 1 + 2^ exp (— (2kx) 2 /2)[1 — (2 kx) 2 } with E(Mq") = \fn/2. 

k=\ 

For / € C[0, oo), let /|[o,i] denotes the restriction of / over [0,1]. Our next 
result is about the weak convergence of the width process Dn ’°^ | [o,i] an d 

cluster process under diffusive scaling. Here and subsequently, 

as is commonly used in statistics, we use the notation X\Y to denote the 
conditional random variable X given Y. 

Theorem 1.3. As n—> oo, we have: 

(i) D { y ) \ m \l {Lm>n} ^V2W+, 

(ii) sup{| pD ( n' 0) (s) - K^’ 0) (s) | : s € [0,1]}1 1 {L(0,0)>n} * 0- 

The following corollary is an immediate consequence of Theorem 1.3. 
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Corollary 1.3.1. Foru> 0, asn^-oo we have: 

(i) y/nP(#C n f 0,0) > y/njou) -A exp(-« 2 /4p 2 ), 

(ii) PEfc=o#ft(0,0) > n 3/2 7 0 u|L(0,0) > n) -A- P(py/2I + > u ). 

Before we proceed to state Theorem 1.4, we recall some results regarding 
random vectors whose distribution functions have regularly varying tails (see 
Resnick [25], page 172). A random vector Z on (0,oo) d with a distribution 
function F has a regularly varying tail if, as n —>■ oo, there exists a sequence 
b n —>• oo such that nP{Z/b n €•}—>• u{-) for some v € M + where M + := {fi : p 
is a nonnegative Radon measure on (0,oo) d }. Here, A denotes vague con¬ 
vergence. It is in this context that Theorem 1.4 obtains a regularly varying 
tail for the distribution of (L(x,t), (#C(x,f)) 2//3 ); which justifies that the 
exponent of Hack’s law is 2/3 for Howard’s model. In addition, we obtain 
a scaling law, with a Hack exponent of 1/2, for the length of the stream, 
vis-a-vis the maximum width of the region of precipitation, that is, 

(5) Anax(0,0) := max{D fc (0,0) : 0 <k< L(0,0)}. 

It should be noted that Leopold and Langbein [18] obtained an exponent of 
0.64 through computer simulations. 


Theorem 1.4. Let E := [0,oo) x [0,oo) \ {(0,0)}. There exist measures 
fi and v on the Borel cr-algebra on E such that for any Borel set B C E we 
have 


( 6 ) 


y/nF 


' (L(0,0),(#C(0,0)) 2 / 3 ) 

n 


€B 




(7) 


y/nP 


(-7/0,0), (D max (0.0)) 1 / 2 ) 


€ B 


n 




with g, and v being given by 



where fj+ and f M + denote the density functions of Iq and Mq" , respectively. 
Moreover, for A,r > 0, we have 


y/nP 


(L(0,0),(#C(0,0)) c 


€ (r, oo) x (A, oo) 


(8) 


n 
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and 

(9) 



, V 7rr 7o ’ 


2 

ot<~, 

2 

1 a> 3 


■y/nP 


(L(0,0),(£> maj ,(0,0)) a ) 

n 


G (r, oo) x (A, oo) 



, \A r 7o' 


1 

ifa< 

1 

lfa > 2 - 


The estimates of the densities fj+ and f M + imply that fi and u are finite 
measures on E. An immediate consequence of the above theorem is the 
following. 


Corollary 1.4.1. As n—> oo for u> 0, we have: 

(i) ^P(#C(0,0) > V2rif'yoPu) -> fo°° t ~ 3/2 ^i+ (ut~ 3/2 ) dt, 

(ii) V^ p (- D max(0,0) > ^/2n7opu) -A ^ 7 =^ f 0 °° t~ 3/2 F M + ('ut“ 1/2 ) dt, 

where Fj+ and F M + are the distribution functions of Iq and Mq , respec¬ 
tively, and Fj+ := 1 — F M + := 1 — ■ 

The proofs of the above theorems are based on a scaling of the process. 
In the next section, we define a dual graph and show that as processes, 
under a suitable scaling, the original and the dual processes converge jointly 
to the Brownian web and its dual in distribution (the double Brownian 
web). This invariance principle is used in Sections 3 and 4 to prove the 
theorems. In this connection, it is worth noting that in Proposition 2.7, 
we have provided an alternate characterization of the dual of Brownian web 
which is of independent interest. This characterization is suitable for proving 
the joint convergence of coalescing noncrossing path family and its dual to 
the double Brownian web and has been used in Theorem 2.9 to achieve the 
required convergence. 

We should mention here that the Brownian web appears as a universal 
scaling limit for various network models (see Fontes et al. [9], Ferrari, Fontes 
and Wu [ 8 ], Coletti, Fontes and Dias [5]). It is reasonable to expect that with 
suitable modifications our method will give similar results in other network 
models. Our results will hold for any network model which admits a dual 
and satisfies (i) conditions listed in Remark 2.1, (ii) the scaled model and its 
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dual converges weakly to the double Brownian web (see Section 2) and (iii) a 
certain sequence of counting random variables are uniformly integrable (see 
Lemma 3.3). In this sense, our result can be considered as a universality 
class result. 

2. Dual process and the double Brownian web. 

2.1. Dual process. For the graph Q, we now describe a dual process such 
that the set of ancestors C(x,t) (defined in the previous section) of a vertex 
(x,t) € V is bounded by two dual paths. The dependency inherent in the 
graph Q implies that, although the cluster is bounded by two dual paths, 
these paths are not given by independent random walks. The dual vertices 
are precisely the mid-points between two consecutive open vertices on each 
horizontal line {y = n}, n € Z with each dual vertex having a unique offspring 
dual vertex in the negative direction of the y-axis. Before giving a formal 
definition, we direct the attention of the reader to Figure 2. 

For ugZ 2 , we define 

:= infjfc : k > 1, (u(l) + fc,u(2)) £7}, 

( 10 ) 

J u := inf {A; : k > 1, (u(l) — fc,u(2)) € V}. 

Next, we define r(u) := (u(l) + J+,u(2)) and l{ u) := (u(l) — J“,u(2)), as 
the first open point to the right (open right neighbour ) and the first open 
point to the left (open left neighbour ) of u, respectively. For (x,t) € V, let 
r(x,t) := (x + J^ xt ^/2,t) and l(x,t ) := (x — J^“^/2,t) denote, respectively, 

the right dual neighbour and the left dual neighbour of (x,t) in the dual 
vertex set. Finally, the dual vertex set is given by 

V := { r(x,t),l(x,t ) : (x,t) £V}. 



Fig. 2. The black points are open vertices, the gray points are the vertices of the dual 
process and the gray (dashed) paths are the dual paths. 
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For a vertex (u,s) € V, let (v,s — 1) € V be such that the straight line 
segment joining ( u,s ) and (v,s— 1) does not cross any edge in Q. The dual 
edges are edges joining all such (u, s ) and (v, s — 1). Formally, for (u, s) € V, 
we define 

a l (u, s ) := sup{z : (z, s — 1) € V, h(z, s — 1)(1) < u}, 

( 11 ) 

a r (it, s ) := inf{z : (z, s — 1) € V, h(z, s — 1)(1) > u} 

and set h(u , s) := (( a l (u , s) + a r (u, s))/2,s — l). Note that ( a r (u , s), s — 1) and 
(■ a\u , s),s— 1) are the nearest vertices in V to the right and left, respectively, 
of the dual vertex h(u, s). Finally, the edge set of the dual graph Q := (V, E) 
is given by 

E := {((u,s),h(u,s)) : ( u,s ) € V}. 


Remark 2.1. Note that the vertex set of the dual graph is a subset of 
x Z. Before we proceed, we list some properties of the graph Q and its 
dual Q. 

( 1 ) Q uniquely specifies the dual graph Q and the dual edges do not 
intersect the original edges. The construction ensures that Q does not contain 
any circuit. 

(2) For (x,t) € V, the cluster C(x,t ) is enclosed within the dual paths 
starting from r(x,t ) and l(x,t). The boundedness of C(x,t) for every (x,t) € 
V implies that these two dual paths coalesce, thus Q is a single tree. 

(3) Since paths starting from any two open vertices in the original graph 
coalesce and the dual edges do not cross the original edges, there is no bi¬ 
infinite path in Q. 

We now obtain a Markov process from the dual paths. Fix (u, s)£f and 
for k > 1, set h k (u,s ) := s)) where h°(u,s) := (u,s). Letting xj. u,s ^ 

denote the first coordinate of h k (u,s ), it may be observed that is a 

function of X^ u,s ^ and the collection of random variables {(RujL^u) : u(2) = 
s — k — 1 G Z}. Thus, by the random mapping representation (see, e.g., Levin, 
Peres, and Wilmer [19]) we have the following. 

Proposition 2.2. For ( u,s ) € V, the process : k > 0} is a time 

homogeneous Markov process. 


Before we proceed, we make the following observations about the tran¬ 
sition probabilities of the Markov process. Let G be a geometric random 
variable taking values in {1,2 ,...}, that is, P(G = l) =p( 1 —py~ l for l > 1. 
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For any u G Z x Z, the random variables and J~ are i.i.d. copies of the 
geometric random variable G independent of -B u . Further, if ui,U 2 € Z 2 are 
such that ui(l) > U2(l) — 1 and ui(2) = U2(2), the random variables 
and Jh, are also independent. Now, for u (ji Z and, v € Z/2, we have 



(u,s) 

1 


>(u,s) 

L 0 


;(v,s) 


U ) F (^-l/2, S -l) ^(u+1/2,s-1) 

= ¥(G 1 -G 2 = 2v), 


2v) 


where G\ and G 2 are i.i.d. copies of G, defined above. If u € Z and v € Z/2, 
we have, using notations from above 


(13) 


P(X 


(«>«) 


- = w|x£ u,a; = u) 


r(lL,s) 

L 0 


= (1 - p)P(Gi - G 2 = 2u) + pP(G = 2u)/2 + pP(G 


—2v)/2, 


where G\ and G 2 are as above. It is therefore obvious that the transition 
probabilities of depend on whether the present state is an integer or 

not. 

From equations (12) and (13), we state the following. 


Proposition 2.3. For any (u, s ) € V, : k > 0} is an L 2 -martingale 

with respect to the filtration J-/ := a({B u , U u : u G Z 2 , u(2) > s — k}). 


2.2. Dual Brownian web. In this section, we briefly describe the dual 
Brownian web W associated with W and present an alternate characteriza¬ 
tion of the dual Brownian web W. 

The Brownian web (studied extensively by Arratia [1, 2], Toth and Werner 
[30], Fontes et al. [9]) may be viewed as a collection of one-dimensional 
coalescing Brownian motions starting from every point in the space time 
plane M 2 . We recall relevant details from Fontes et al. [9]. 

Let M 2 denote the completion of the space time plane M 2 with respect to 
the metric 


p((xi,ti),(x 2 ,t 2 )) ■■= |tanh(ti) — tanh(t 2 )| V 


tanh(xi) 

1 + |ti| 


tanh(x 2 ) 

i + M 


As a topological space M 2 can be identified with the continuous image of 
[—oo, oo] 2 under a map that identifies the line [—oo, oo] x {oo} with the point 
(*, oo), and the line [—oo, oo] x {—oo} with the point (*, —oo). A path n in M 2 
with starting time a 7r € [— 00 , 00 ] is a mapping it : [ov,oo] —» [— 00 , 00 ] U {*} 
such that 7r(oo) = * and, when ov = — 00 , 7r(— 00 ) = *. Also t i-A (n(t),t) is a 
continuous map from [ov,oo] to (M 2 ,p). We then define II to be the space 
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of all paths in M[? with all possible starting times in [— 00 , 00 ]. The following 
metric, for € II 


(in^DT 2 ) := max-] |tanh(fj 7ri ) — tanh(cr 7r2 )|, 

tanh(7Ti (t V ovj) tanh(7T2(f V cr^)) 


sup 

Act7t 2 


1 + 


1 + 


makes II a complete, separable metric space. 


Remark 2.4. Convergence in this metric can be described as locally 
uniform convergence of paths as well as convergence of starting times. There¬ 
fore, for any e > 0 and m > 0, we can choose £i(= /(e,m)) > 0 such that 
for 7ri,7T2 € II with {(7Tj(i),t) : t € [cr^, m ]} C [— m,m\ x [— m,m\ for i = 
1 , 2 , dn( 7 Ti, 7 r 2 ) < £1 implies that ||( 71 - 1 ( 0 -^), cr^) - (7T 2 (cr^- 2 ), cr^- 2 )|| 2 <£ and 
sup{17r 1 (f) — 7T2(f)| : t € [maxjoVj, av 2 },77i]} < £. We will use this later several 
times. 


Let 7L be the space of compact subsets of (II, du) equipped with the 
Hausdorff metric d The Brownian web W is a random variable taking 
values in the complete separable metric space (TL, d%). 

Before introducing the dual Brownian web, we require a similar metric 
space on the collection of backward paths. As in the definition of II, let 
II be the collection of all paths 7r with starting time cr# € [— 00 , 00 ] such 
that 7T : [— 00 ,a*] -A [— 00 , 00 ] U {*} with 7r(— 00 ) = * and, when a # = + 00 , 
7r(oo) = *. As earlier t i-A (7(t),t) is a continuous map from [— 00 ,( 7 #] to 
(Me, p) ■ We equip II with the metric 


dyr (iti,7T2) := max< |tanh(cr# 1 ) — tanh(cr# 2 )|, 




sup 


tanh(7Ti(t A cr^)) tanh(7T2(f A cr# 2 )) 


1 + 


1 + 


making (II, dg) a complete, separable metric space. The complete separable 
metric space of compact sets of paths of II is denoted by (71, d-q), where dg 
is the Hausdorff metric on TL, and let B ^ be the corresponding Borel a held. 


2.3. Properties of (W,W). The Brownian web and its dual (W,W) is 
a (HxH,B h x B^j) valued random variable such that W and W uniquely 

determine each other almost surely with W being equally distributed as —W, 
the Brownian web rotated 180° about the origin. The interaction between 
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the paths in W and W is that of Skorohod reflection (see Soucaliuc, Toth 
and Werner [28]). 

We introduce some notation to study the sets {7r(t + s) : tt G W, a n < t} 
and {tt (t — s) : tt € W, cr% > t}. For a (H, By) valued random variable K and 
t € R, let K t_ := {n : n € K and a n < t}. Similarly, for a (' H,B valued 
random variable K and let K t+ := {it: tt £ K and > t}. For ti,t 2 € 

R, t 2 > t\ and a (H, By) valued random variable K, define 


(14) 


:={7r(f 2 ) '.TreK* 1 ,n(t 2 ) € [0,1]}; 
^K(h,t 2 ) := #A4^(ti,t 2 ), 


that is, ^(ti, £ 2 ) denotes the number of distinct points in [0,1] x f 2 which are 
on some path in K fl ~ . We note that for t > 0, M.\\>(to, to + t) = Ayv(to- 1; 0,1) 
as defined in Sun and Swart [29] . It is known that for all t > 0 the random 
variable +t) is finite almost surely (see (£}) in Theorem 1.3 in Sun 

and Swart [29]) with 

(15) E(£ w (t 0 ,t 0 + t)) = -^=. 

y/TTt 

Moreover, from the known properties of (W\W) the proof of the following 
proposition is straightforward (for details, see Roy, Saha and Sarkar [27]). 


Proposition 2.5. For any to < t \ , almost surely we have: 

(i) My\>(to,ti) n Q = 0 ; 

(ii) each point in My\>{to,t{) is of type (1,1); 

(iii) for each x € A4y\>(to,ti), there exists 7Ti,7t 2 € W with a ni < to, a^ 2 > 
to and 7Ti(ti) = 7r 2 (ti) = x; 

(iv) for each x € M\\i(to,ti), there exist exactly two paths Tti x,tl ^ and 
Tt( x,tl ) in W starting from (x,t\) with itl- x,tl \t) > 7 Tj x,tl \t) for all [to? ^ 1 )- 

There are several ways to construct W from W. In this paper, we follow 
the wedge characterization provided by Sun and Swart [29]. For 7r r ,7r z € W 
with coalescing time t nT,n and 7r r (max{ a n r, cr n i}) > tt 1 (max.{a n r, a n i}), the 
wedge with right boundary tt t and left boundary tt 1 , is an open set in R 2 
given by 

A = A(TT r ,TT l ) 

( 16 ) 

:= {(y,s) : ma x{a n i,a 7r r} <s< t^'* ,n l (s) <y< 7r r (s)}. 

A path tt £ II, is said to enter the wedge A from outside if there exist t\ and 
t 2 with o^>ti> t 2 such that (7r(ti), ti) ^ A and (7r(f 2 ),t 2 ) € A. 
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From Theorem 1.9 in Sun and Swart [29] , it follows that the dual Brownian 
web VV associated with the Brownian web VV satisfies the following wedge 
characterization. 

Theorem 2.6. Let (W ; ,VV) be a Brownian web and its dual. Then al¬ 
most surely 

VV = {fr: 7r € n and does not enter any wedge in VV from outside}. 

Because of Theorem 2.6, for a (TL x U,B n x By) valued random variable 

(W,Z) to show that Z = VV, it suffices to check that Z satisfies the wedge 
condition. Here we present an alternate condition which is easier to check. 

Proposition 2.7. Let (W,Z) be a (TL x TL,B% x B^) valued random 
variable such that: 

(1) for any deterministic (x,t) G M 2 , there exists a path n( x,t ^ G Z starting 
at (x,t) and going backward in time almost surely; 

(2) paths in Z do not cross paths in VV almost surely, that is, there 
does not exist any tt G VV, iG2 and t\,t 2 G (a n ,a^) such that (vf(fi) — 
7r(ti))(7r(t2) — 7r(t2)) < 0 almost surely; 

(3) paths in Z and paths in VV do not coincide over any time interval 
almost surely, that is, for any it G VV and tt G Z and for no pair of points 
t\ < t 2 with <7,,- < t\ < t 2 < Utt we have tt( t) = tt( t) for all t G [H,^] almost 
surely. 

Then Z = VV almost surely. 

Proof. From conditions (2) and (3), we have that tt G Z does not enter 
any wedge in VV from outside. Hence, Z C VV. The argument for VV C Z 
follows from the fish-trap technique introduced in the proof of Lemma 4.7 

of Sun and Swart [29]. It shows that VV C Z almost surely for any (TL,Bu) 
valued random variable Z satisfying (i) paths in Z do not cross paths is VV 
and (ii) for any deterministic countable dense set, there exist paths in Z 
starting from every point of that dense set (for details, see Roy, Saha and 
Sarkar [27]). □ 

2.4. Convergence to the double Brownian web. For any (x,t) G V, the 
path 7 r( x P in the random graph Q is obtained as the piecewise linear function 
7 r( x,t ) : [t,oo) —>■ M with 7 r( x,t \t + k) = h k (x,t)( 1) for every k > 0 and 7 
being linear in the interval [t + k,t + k + 1], Similarly, for (x,t) G V, the 
dual path 7 r^’*) is the piecewise linear function 7r( x,t ) : (— 00 , t] — > M with 
TT^ x,t \t — k) = h k (x,t)( 1) for every k > 0 and Tt^ x,t ^ being linear in the interval 
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[t — k — l,t — k\. Let X := {n&F) : (x, t) € V} and X := {7 :(i,t)gV} be 
the collection of all possible paths and dual paths admitted by Q and Q. 

For a given 7 > 0 and a path 7 r with starting time 07 -, the scaled path 
717 ( 7 ) : [ 07 -/ 71 , 00 ] -A [— 00 , 00 ] is given by 717 ( 7 )^) = itint) / (y/n'y) for each 
n > 1. Thus, the starting time of the scaled path 717 ( 7 ) is 07 n ( 7) = 07 /n. Sim¬ 
ilarly, for the backward path 7 r, the scaled version is 717 ,( 7 ) : [— 00 , 07 -/ 71 ] -A 
[— 00 , 00 ] given by 7 r n ( 7 )(t) = 7 r(nt)/ {y/n'y) for each n > 1. For each n > 1, let 
X n = X n {^) := { 7 T^’ <) ( 7 ) : ( x,t ) € V} and X n = X n ( 7 ) := { 7 r, ( /’ <) ( 7 ) : (x,f) € 
V} be the collections of all the nth order diffusively scaled paths and dual 
paths, respectively. 

The closure X n {^/) of X n {j) in (II,dn) and the closure X n {p/) of X n {p /) 
in (II,dg) are (LL&h) and (TL,By) valued random variables, respectively. 
Coletti, Fontes and Dias [5] showed the following. 

Theorem 2.8. For 70 := 7 o (p) as in Theorem 1.2, as n->- 00 , X n ( 70 ) 
converges weakly to the standard Brownian web W. 

Our main result is the joint invariance principle for {(T n ( 7 o),T n ( 7 o)) : 
n > 1 } considered as (TL x 'H,B-n x By) valued random variables. 

Theorem 2.9. {(T n ( 7 o),T n ( 7 o)) : n > 0} converges weakly to (W,W) 

as n —> 00 . 

We require the following propositions to prove Theorem 2.9. We say that 
( WOT (n) : u < t} is a Brownian motion going back in time if (t — s ) := 

W(t + s), s > 0 where {W (it) : u > t} is a Brownian motion with W ( t ) = x. 

Proposition 2.10. For any deterministic point ( x,t ) gM 2 , there exists 
a sequence of paths 9n~ t ' 1 G X n { 70 ) which converges in distribution to W^ x,t \ 

Proof. For any ( x,t ) G M 2 fix t n = [nt\ and x n = maxjL-y/nyoxJ + j : 
j < 0, ([yTiyoxJ + j,t n ) G V}. Let G X n ( 70 ) be the scaling of the path 

TT^U) G % 

Since Q is invariant under translation by lattice points and Q is uniquely 
determined by Q , the conditional distribution of {( 27 , t n ) + h ?( 0 , 0 ) : j > 0 } 
given (0,0) G V is the same as that of {hi ( 27 , t n ) : j > 0}. We observe that 
(x n /(y/n^o),t n /n) —> ( x,t ) as n A 00 almost surely. Hence, it suffices to 
prove that the scaled dual path starting from (0,0) given (0,0) € V converges 
in distribution to W^ 0,0 ). 
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From Proposition 2.3, we see that xj°’ 0) = hP(0, 0)(1) is an L 2 martingale 
with respect to the filtration o ({-£>(.,<,), Ut z ^ : z £Z,s> —k}). Let 

%(«) : = + (-Xj'J? - Xj m )(uS 2 n - s))/(s 2 j+1 - S 2 )] 

for u € [ 0 , oo) and s 2 < us 2 < s 2 +1 , where s 2 = J2j=i ^((xj 0 ’ 0 ^ - xj 0 ’^) 2 ). 
We know r/ n converges in distribution to a standard Brownian motion (see 
Theorem 3, [4]). Since Sn/( n 7o) ~^ 1> it can be seen that sup, ug [ 0 ^ \r] n (u ) — 

u)| —>■ 0 in probability for any M > 0. So by Slutsky’s theorem, we 

conclude that 6n’°^ converges in distribution to a standard Brownian motion 
going backward in time. □ 

The next result helps in estimating the probability that a direct path and 
a dual path stay close to each other for some time period. Given m € N and 
e, 8 > 0 , we define the event 

Bn = B £ n (5,m ) 

:= {there exist 7 r™, tt% , £ X n such that 077 , 0 V" < 0, 

07 -n < |_w< 5 J/n, 7 r”( 0 ) € [—m,m], | 7 r”( 0 ) — ( 0 )| < s, with 

ni (\n5 \ / n) tt 2 ([nd \ / n) and \^i{[n5\/n) — ir^ (|_n<5J/n)| < e, with 

<(2 |n<5J /n) / vr^ (2 [nS\ /n )}. 

Lemma 2.11. For am/ m € N and e,5 > 0, we hare 

P(F^(h, m)) < Ci(S,m)e, 

where Ci(S,m) is a positive constant, depending only on 5 and m. 

Proof. Let D e n be the unsealed version of the event B e n , that is, 

D £ n := {there exist (x,0), (y, 0), (z, |_n<5J) £ V such that 

x £ [-m^n'y 0 ,m\/n'yo], \x — y\ < \fne 7o and (x, 0 ) / (y, 0 ), 

\h^ (®,0)(1) -z\< v^e 7 o ,h 2 ^ (x,0) + (z, Ln<5J)}. 

On the event D £ n there exists l £ [—m^/n'yo, m^/n r yo] HZ such that the un¬ 
sealed paths starting from (l, 0) and (l + 1,0) (as in Figure 3) do not meet in 
time [n5\ —an event which occurs with probability at most C 2 /VnS for some 
constant C 2 > 0 (see Theorem 4 of Coletti, Fontes and Dias [5]). Supposing 
h,L n<5 J (Z, 0 )( 1 ) = k, two unsealed paths, one starting from a vertex [^/ne 70 J 
distance to the left of k and the other starting from a vertex |_\/™ e 7 oJ dis¬ 
tance to the right of k, do not meet in time has a probability at most 
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Fig. 3. The vertices (£,0) and (Z H-1,0) and the corresponding vertex ( k, |_n(5J) as required 
in the proof of Lemma 2.11. 


C 2 2y / n£7o/v / n<5 for all feEZ. Thus, summing over all possibilities of l and 
k and using Markov property we have 


( 2m v /n70 

u 1> M (' ,0)(1) = k^h lnSi {l+ 1,0){1) and 

l=—2m^fn r )o AiGZ 

h ^ nS J (k - L\/n£7oJ, |n<5J) / h^ nSi (k + |_\/™£7oJ, Ln5J)} ) 


2m v /n7o r 

< E 2 )v° E(*> w 1 )=* hln5i V +0)(!)> 

l=—2my/n r yo ^ fceZ 


2my/n 10 

< E ^^P £ IP{^ Ln ' 5J (^0)(l)//i Ln ' 5J (Z + l,0)(l)} 

l=—2my/n r )Q 
2m^n7o 

£ £ 

l=—2myfn')o 

< Ci(<5, m)e. 


Vn5 

2C < 2 y / n£7o C 2 


□ 










HACK’S LAW IN A DRAINAGE NETWORK MODEL 


17 


Proof of Theorem 2.9. Since X consists of noncrossing paths only, 
Proposition 2.10 implies the tightness of the family {X n : n > 1} (see Propo¬ 
sition B.2 in the Appendix of Fontes et al. [9]). The joint family {(X n ,X n ) : 
n > 1} is tight since each of the two marginal families is tight. To prove 
Theorem 2.9, it suffices to show that for any subsequential limit (W,i?) of 

{(X n , X n ) : n > 1}, the random variable Z satisfies the conditions given in 
Proposition 2.7. 

Consider a convergent subsequence of {(X n , X n ) : n > 1} such that (W, Z) 
is its weak limit and by Skorohod’s representation theorem, we may assume 
that the convergence happens almost surely. For ease of notation, we denote 
the convergent subsequence by itself. 

From Proposition 2.10, it follows that for any deterministic (x,f) € M 2 
there exists a path ir £ Z starting at ( x,t ) going backward in time almost 
surely. 

Since (X n ,X n ) converges to ( W,Z ) almost surely, if a dual path in Z 
crosses a path in W, there exists a dual path in X n which crosses a path in 
X n , for some n > 1, yielding a contradiction. Hence, the paths in Z do not 
cross paths in W almost surely (for details, see Roy, Saha and Sarkar [27]). 

Now, to prove that condition (3) in Proposition 2.7 is satisfied, we define 
the following event: for 5 > 0 and positive integer m > 1, let 

A(6,m ) := {there exist paths it G W and tt € Z with € (— m,m), 

and there exists to such that cr n < to < to + 6 < a ft, 

and — m < ir(t) = 7r(f) < m for all t € [to, to + £]}■ 

It is enough to show that for any fixed 5 > 0 and for m > 1, we have 
¥{A(d,m)) =0. 

We present here the idea of the proof; more details are available in Roy, 
Saha and Sarkar [27]. Fix e > 0. Since we are in a setup where the scaled 
paths converge almost surely, for all large n there exist ir " € X n and n n € X n 
within e distance of ir and n, respectively. Using the fact that a dual vertex 
lies in the middle of two open vertices and the forward paths cannot cross 
the dual paths, it follows that for all large n there exist vr^vrg € X n such 
that: 

(a) max{|7rf (a n %) - tt£(ct^)!, |vrf (a^) - vif (cr^)!) < 4e; 

(b) n^(a n n + 6/3) / ^(ovj + 6/3) and + 6/3) / ^(cr^n + 6/3). 

This gives us that A(6, m) C lirriinf' rwoo ^ B//{6/3 1 2m]j). 

Here, B^ £ (6/3,2m;j) is a translation of the event B^ £ (6/3,2m), consid¬ 
ered in Lemma 2.11; translated such that the starting time of the paths tt \ 
and 7r 2 are shifted by — m + j[n6/3\/n (see Figure 4). 
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S 



Fig. 4. The event A(5, m). The bold paths are from (W, W) and the approximating dashed 
paths are from (X n ,X n ). 


By translation invariance of our model and Lemma 2.11, for all n > 1 we 
have P(L>^ £ ((5/3,2m; j)) < 4Ci(<5/3, 2m)e. This completes the proof. □ 

3. Proof of Theorem 1.2. Let £ := £w(0,1) and £ n := £^(0,1) be as 
defined in (14). The proof of Theorem 1.2 follows from the following propo¬ 
sition. 

Proposition 3.1. E[£ n ] -aE[£] as n^oo. 

We first complete the proof of Theorem 1.2 assuming Proposition 3.1. 


Proof of Theorem 1.2. Using the translation invariance of our model, 
we have 


LvWoJ 


v^7o 


v / n7 0 P(L(0,0) > n) = ^ E(l {i ( M)>n }) 


Lv^ToJ +1 

_ 1 


k =0 



->E(£) = -= 

VT 


as no oo. 


This proves Theorem 1.2. □ 


Proposition 3.1 will be proved through a sequence of lemmas. 
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To state the next lemma, we recall from Theorem 2.9 that (X n ,X n ) =>- 
(W, W) as n —>■ oo. Using Skorohod’s representation theorem, we assume 

that we are working on a probability space where d Hx ^((X n , X n ), (VV, W)) —>• 
0 almost surely as n —>• oo. 

Lemma 3.2. For ti > to; we hare 

^(£a„(^o>U) ^€w(to,h) for infinitely many n) = 0 . 

Proof. We prove the lemma for to = 0 and t\ = 1, that is, for £ n = 
(0,1) aR d £yy(0, 1), the proof for general to,ti being similar. First, we 
show that, for all k > 0, 

(17) liininf l { ^> fc} > l {? > fc} almost surely. 

Indeed, for k = 0, both l{£ n >fc} and l{£>fc} equal 1. For k > 1, (17) follows 

from almost sure convergence of (X n , X n ) to (W, W) and from the properties 
of the set Ad>v(0,1) as described in Proposition 2.5. 

To complete the proof, we need to show that P(limsup n _ ) . 00 {£ n > £}) = 0. 
This is equivalent to showing that P(flg) = 0 for all k> 0, where 

fig := {u : £„(cj) > £(w) = k for infinitely many re}. 

Consider k = 0 first. From Proposition 2.5, it follows that on the event £ = 0, 
almost surely we can obtain 7 := 7 ( 0 ;) > 0 such that .M>v( 0 , 1 ) n (— 7,1 + 

7 ) = 0. From the almost sure convergence of (X n ,X n ) to (W, W), we have 

m° 0 )=o. 

For k > 0, on the event fig we show a forward path n € W coincides with 
a dual path it € W for a positive time which leads to a contradiction. From 
Proposition 2.5, it follows that given 7 > 0, there exist mo € N and so € 
(l/m 0 ,l) such that P(£ w (l/m 0 ,1) = £w(l/m 0 , s 0 ) = £w( 0 ,l) = k) > 1 - r], 
that is, the paths leading to any single point considered in Adyv( 0 , 1 ) = 
X4 yy(l/77io, 1 ) have coalesced before time so- Fi x 0 < e < 1 /mo such that 
(x — s, x + e) C (0,1) for all x € 1) and the e-tubes around the k 

paths contributing to Adw;(so,l), viz., 717(t),..., 7Tk{t), t € [so,l], given by 

T\ := {(x. t ) : 7 Ti(t) — e < x < 1 r,(t) + e, so < t < 1 } for i = 1 ,..., k, 

are disjoint. Since we have almost sure convergence on the event fig, there ex¬ 
ists no such that one of the k tubes must contain at least two paths, 7r”°, irlf 0 
(say) of X no which do not coalesce by time 1 . From the construction of dual 

paths, it follows that there exists at least one dual path k n ° G X nQ lying be¬ 
tween 7 T™° and 7 T 2 ° for t € [sq, 1 ], and hence we must have an approximating 
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7 f E W 1+ close to 7r n ° for t E [so, 1]. Since we have only finitely many disjoint 
k tubes, taking e —>• 0 and using compactness of W we obtain that there 
exists 7r G W such that 7r(t) = 7Tj(f) for t E [so, 1] and for some 1 < z < k. 
This violates the property of Brownian web and its dual that they do not 
spend positive Lebesgue time together. Hence, P(fi§) = 0 for all k> 0 and 
this completes the proof of the lemma. □ 

Lemma 3.2 immediately gives the following corollary. 

Corollary 3.2.1. As n—>■ oo, £ n converges in distribution to £. 


Corollary 3.2.1 along with the following lemma completes the proof of 
Proposition 3.1. 

Lemma 3.3. The family {£ n : n E N} is uniformly integrable. 


Proof. For me N, let 

A ra =[-m,m] 2 nZ 2 and := {(0,1), (0,-1), (1,1), (1,-1)} A ’ 


We assign the product probability measure P 7 whose marginals for u e K m 
are given by 


P'{C:C(u) = (a,6)} 


P 

2 ’ 


(!~P) 

2 


for a = 1 and b E {1, — 1}, 
for a = 0 and b E {1, — 1}. 


P 7 is the measure induced by the random variables {(H u , U u ) : u E K m }. 

For £ E and for I\ C K m , the K cylinder of £ is given by C{Q,I\) := 
{£ 7 : £ 7 ( u ) = £( u ) for all u E K }. For any two events A,B C Q m , let 


A\AB := {£ : there exists K = K{(f) C K m such that C(£, K) C A, 
and C((, K') C B for K' = K m \ K} 


denote the disjoint occurrence of A and B. Note that this definition is asso¬ 
ciative, that is, for any A,B,C C we have (A\AB)\3C = A\A(B\I\C). 

Let 


F™ := {there exist (u±,n), (u 2 ,n) E V with 0 < u\ < U 2 < y/n^o and 
(v[,l), (y l 2 ,1) £ V for all 0 < l < n such that 
—m < v\ < h l (ui,n)(l) < h l {u 2 ,n){\) <v l 2 < m}, 


E™(k) := {for 1 < i / j < k, there exists (xj,0) E V with 

h n (xi, 0)(1) E [0, \/n r fo\,h n (xi, 0) / h n (xj, 0), h l (xi, 0)(1) E [-m,m\ 
for all 0 < l < n}. 
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(ui,n) (u 2 ,n) (u, 3 ,n) ( 114 , n) (u 5 ,n) 



Fig. 5. The event E™ (6). 


We claim that for all k> 2, 

(18) E™(3k) C F™UF™U ■ • • OF? . 

V v y 

k times 

We prove it for k = 2. For general k, the proof is similar. Let (ui,n) € 
V, 1 < i < 5 and (xj, 0) € V, 1 < % < 6 be as in Figure 5. The region explored 
to obtain the vertex h^(ui,n) for 1 < j < n is contained in (J)]l 0 1 [h*(xi, 0 )(l), 
h l (xi+\, 0)(1)] x {Z}. Thus, the regions explored to obtain the dual paths 
starting from (ui,n), (u 2 ,n) and the dual paths starting from ( 7 / 4 , 71 ), {u§,n) 
are disjoint (see Figure 5). Hence, it follows that _E™( 6 ) C F™OF™. 

Since the event E™(k) is monotonic in m, from (18) we get 

P(£ n > 3k) = pf lim F™(3k))= lim P(F™(3k)) 

\ra—> oo / ra—>-oo 

< lim P(F™D • • • OF™) = lim P '{F™U-■-UF™). 

m—>• 00 m—>-oo 

Applying the BKR inequality (see Reimer [24]), we get 

(19) P(£„ > 3k) < lim (P'(F n m )) fc = (p( lim F™)Y = (P (F n )) k , 

m—>•oo V \m —>oo / / 

where F n := {there exist (ui,n), ( 7 / 2 , 71 ) G V with 0 < u\ < 112 < y^To such 
that h n (ui,n) / (7/2,71)}. 

For any (x,t) G M 2 fix t n = |_rafj and x n = max{ \_s/vrfyx\ + j : j < 0, 
(LVnyoxJ +j, i n ) G V}. Let G A n ( 7 o) be the scaling of the path 7 T^ n,tn ) G 
A. Dehne 

F' n := {0) ) 0 ’ 1 ^ and do not coalesce in time 1}. 
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We observe that F n <Z F' n . Now P(i^) converges to the probability that 
two independent Brownian motions starting at a distance 1 from each other 
do not meet by time 1. Since lim^oo P(F^) < 1, the family j£ n : n € N} is 
uniformly integrable. □ 

Remark 3.4. It is to be noted that Newman, Ravishankar and Sun [23] 
also used ideas of negative correlation to establish the weak convergence of 
M-x as a point process on M for a more general setup where paths can cross 
each other. In our case, the negative correlation ideas come in a much less 
essential manner only to establish uniform integrability as the noncrossing 
nature of paths enable us to obtain Corollary 3.2.1. 


4. Proofs of Theorems 1.3 and 1.4. In this section, we prove Theorems 
1.3 and 1.4. The main idea of the proof is that the horizontal distance 
between the dual paths n f ( x ^ and (see Figure 6) form a Brownian 

excursion process after scaling. The cluster C(x,t ) being enclosed between 
these two paths, its size is related to the area under the Brownian excursion. 

For a formal proof, we need to introduce some notation. For r > 0, let 
S T , S T+ : C[0, oo) ->R be defined by S T (f ) := inf jt > 0 : f(t + s) > f(t ) for 
all 0 < s < r} and S T+ (/) := infjf > 0 : f(t + s) > /(f) for all 0 < s < r}. Let 
T T+ : CjO, oo) —>• C[0, oo) be the map given by 


(20) T T+ (/)(«) 


f(S x+ +s)-f(S x+ ), if S T+ < oo, 

f(s), otherwise. 



Fig. 6. The two dual paths and enclose the cluster C(x,t). These dual 

paths after scaling are each Brownian paths. 
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For a Brownian motion W with W{ 0) = 0, we define W T = T T+ (W). From 
Bolthausen [3] , we have S T+ = S T < oo almost surely under the measure in¬ 
duced by W on C[0,oo) and PF 1 |[ 0 ,i] = W + where W + is the standard Brow¬ 
nian meander process defined in (3). From the scaling property of Brownian 

motion, it follows that { W T (s) :sE[0,r]} = {\/tW + (s/t) :sG[0,t]}. Dur- 
rett, Iglehart and Miller [7] (Theorem 2.1) proved that W|l{ min ^ w(s)>-e} 
W + as e 10. Using this result and the scaling property of W T , given above, 
straightforward calculations imply the following lemma and its corollary (for 
details, see Roy, Saha and Sarkar [27]). 

Lemma 4.1. For r > 0 considering W as a standard Brownian motion 
on [0, oo) starting from 0, we have W |l{min te[0 r , w(t)>-i/n} => W T as n —>■ oo. 


Define W T as the process on C[0, oo) given by 


W T (t) 


W T (t), _ if0<t<r, 

W t {t) + W(t — t), otherwise, 


where W is a Brownian motion on [0,oo), independent of W T , with W(0) = 
0. For / € C[0, oo), let tf := inf{s > 0 : f(s) = 0} with tf = oo if f(s) / 0 for 
all s > 0. Consider the mapping H : C[0,oo) —» ^[0,00) given by H(f)(t) : = 
1 {t<tf}f(t)- We define W +,T = H(W T ). A similar argument as that of Lemma 
4.1 gives us the following corollary. 


Corollary 4.1.1. Forr> 0, we have, W T = W T andW + ’ T = H(W T ). 


Let A C C[0, 00 ) be such that 

A := {/ € C[0, 00 ) '.tf < 00 and for every e > 0 there exists 

( 21 ) 

s € + e) with f(s) < 0}. 

From Corollary 4.1.1, it follows that P(T / F T € A) = 1. Hence, H is continuous 
almost surely under the measure induced by W T on (7[0,oo). 

Next, we obtain the distribution of / 0 °° W +,T (t) dt. 

Lemma 4.2. For r, A > 0, we have 

p(jf W + ' T {t)dt >\j = ^ t-^F^iXt-^dt. 

Proof. We give here a straightforward proof using random walk. Let 
{5 n : n > 0} be a symmetric random walk with variance 1 starting at So = 0. 
Since P (W T € A) = 1, minor modification of the argument used to prove 
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Lemmas 2.4 and 2.5, Bolthausen [3] shows that H oT T is almost surely 
continuous under the measure induced by W on C[0,oo) (for details, see 
Roy, Saha and Sarkar [27]). From Donsker’s invariance principle and from 
the continuous mapping theorem, it follows that for A > 0, a continuity point 
of / 0 °° W +,T (t) dt, we have 

F (/ W+ ’ T ^ dt> A ) = n^c ¥ (J 0 H ( TT+ ( Y nW)dt>^j, 

where 

(22) y„(t):=^ + M_^( Si+1 - St ) for l<t<l±I 

\ n \ n n n 


A similar argument as in Lemma 3.1 of Bolthausen [3] gives us that (for 
details, see Roy, Saha and Sarkar [27]) 


</: 


H(T T+ (Y n )){t)dt> A 


/ H(Y n ){t)dt>\ min Y n [t) > 0,to > nT 

1 0 £6[0,t] 


where to := inf{n > 0 : S n = 0} is the first return time to 0 of the random 
walk. Hence for A > 0, a continuity point of W +,T , we obtain 


f 

Jo 


W + ’ T (t)dt> A 


= lim P 

n—>■ oo 


/ H ( Y n )(t ) dt > A min Y n (t) > 0, to > nT 

1 0 i£[0,r] 


lim \ 

n—>oo * 


n 3 / 2 P(io = [nrj + j) 


n(y/n¥(to> nr)) 


x P 


r oo 

/ H(Y n )(t) 
Jo 


dt>\ 


min Y n {t ) > 0,t 0 = [nrj + j 
te[o ,t] 


= lim 


9n(t)fn(t)dt, 


n^oo y/nP(t 0 > nr) 7 L nrJ/n 

where for f > [nT\/n, f n (t ) = P(/ 0 °° H(Y n )(u) du > A] min te [ 0ir ] Y n (t) > 0,t 0 = 
|_ntj + 1) and g n (t ) = n 3 / 2 P(to = [ntj + 1). It is known that (see Kaigh [16]) 


lim \/nP(to > n) 

n—t oo 



lim n 3,/2 P(to = n) 

n^-oo 



Hence, from Theorem 2.6 Kaigh [16] together with the continuous map¬ 
ping theorem and the scaling property of the Brownian motion we have 
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P(/ 0 °° W +,T (t ) dt > A) = ^ / r °° t ^F I+ (A t 3 / 2 ) dt. Finally, Iq being a con¬ 
tinuous random variable (see Louchard and Janson [20]), it follows that the 
random variable J 0 °° W +,T (t ) dt is continuous. This completes the proof. □ 


4.1. Proof of Theorem 1.3. Recall that r(x,t ) and l(x,t ) denote the 
right and left dual neighbours, respectively, of (x,t) G V. Let Df.(x,t) := 
h k (r(x,t))( 1) — h k (l(x,t))( 1) where h is as defined after (11). Consider the 
continuous function G C[0,oo) given by 


(23) 


Dt%) := - 4(x,t» 

7ov n 7oV n 


. k k + 1 

for — < s <-. 

n n 


Fix r > 0. For an TL x TL valued random variable ( I \, I \) and for x € 
A4^(0,r) let be defined as 


_ f 7r, if <Tjr = r and there is no fri G iC r+ with x < 7tt(t) < 7r(r), 

\ 7Tq, otherwise, 

where 7To denotes the constant zero function with a^ 0 = r. In other words, 
7Tr l ’ r ) G iT r+ is such that among all 7r G iL r+ , is closest to (x,r) on 

the right. Similarly, fry ,r is defined as the path closest to (x,r) on the left. 

For 7r G II with > r, let o(7r) G C[0, oo) be given by q(n)(t) := 7r(r — t) 
for t > 0. Fix / G C b [ 0, oo) and define 

K (A,A)( T ’/) : = X] f(9^i x ’ r) )-9(^ X ’ t) ))- 

xGMk(.0,t) 


Let k(t,/) := K (vy ^ } (r, /), and K n (r,f) := n- ~ ,(t,/). Comparing with 




the definitions introduced in (14), for m/ = sup{|/(s)| : s G [0,oo)} we have 

(24) k(tJ) <m/^w(0,r),/c n (r,/) <m/^ n (0,r) foralln>l. 

From Proposition 2.5, we know that for each x G A4>v(0, r), there exist 
7Tr Z ’ r \ G W both starting from (x,r) with ^’^(O) > 7T^’ t ^(0). 

The following lemma is the main tool for establishing Theorem 1.3 and 
Theorem 1.4. 


Lemma 4.3. For r > 0 and f G Cf,[0, oo), we have 
(25) lim E[/c n (r, /)] =E[k(t,/)]. 
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Proof. From (24) and Lemma 3.3, it follows that the family {n n (r, /) : 
n € N} is uniformly integrable. Hence, it suffices to show that At n (r,/) con¬ 
verges in distribution to At (r, /) as n —> oo. We assume that we are working 


on a probability space such that (X n , X n ) converges to (W, W) almost surely 
in (Ti x From Lemma 3.2, we have linp^oo^ (0, r) = 6v(0,t) 

almost surely, and hence from (24) for ^w(0,r) = 0, we have At n (r,/) = 
At (r, /) = 0 for all n large. Next, we consider the case £w(0, r) = k > 1. Sup¬ 
pose A4w(0,r) = {xi,... ,Xk}. From Lemma 3.2, we have that (0, r) = 


{x^,..., x^} for all large n and lim^^oo x™ = X, for all 1 <i <k. Fix T > 0. 
To complete the proof, it is enough to show that sup{|7fr X,,T ^(r — s) — 

7r^ 1 ,t \t — s)| V |7r^’ r ^(r — s) — ’"^(t — s)| : s € [0, r + T]} —>• 0 as n —> oo 

for all 1 < i < k. 

We observe that for yi € there exists 77^’°) € W 

such that 7r^ i,0 )(r) = Xj. We choose s = s(oj) > 0 so that for all 1 < i < k: 


(a) (xi — e, Xj + e) C (0,1), (xj — 2s, Xj + 2e) n A4w(0, r) = {x^} and 

(b) (# i,T) (0) - 7T^’°)(0)) A (7r^’°)(0) - 7T; (:Ei ’ r) (0)) > 2s. 

Let no = no(w) be such that, for all n>riQ\ 

(i) 6 f„( 0 ,t)=£ w ( 0 ,t) and 

(ii) for all 1 <i <k there exist fry”, 77 2,n € and 77™ € such that 

sup{|7r J 1,n (r - s) — - s)| V |7r J 2,n (r - s) — - s)| V |7r"(r — s) — 

-jrivifi )(r — s)| : s € [0,r + T]} < e. 

The choice of no ensures that (0, r) n (xj — e, Xi + e) = {x™}. Since there 
exist only two dual paths starting from (xj,r), because of the uniqueness of 
x™ in the interval (xj — e, Xj + e) and the noncrossing nature of our paths we 

must have 77 ,T \r — s) = 77 *’" (r — s) and 7r| x ’ ’^(t — s) = 7r 2,n (r — s) for all 
s € [0, r + T] and for all n > n q (for details, see Roy, Saha and Sarkar [27]). 
Since T > 0 is chosen arbitrarily, this completes the proof. □ 


The next lemma calculates E[av(t,/)]. 

Lemma 4.4. For r > 0 and f € C&[0, oo), we have 

EKt,/)]=E(/(V2W+’"))MF. 

Proof. Let I n C {0,1,... ,n — 1} given by I n := {i: 0 < i < n— 1, ^ l l n ^ T \ 
7r((* +1 )/ n,r ) e W such that 7r(* //n,r ) (0) < 7r^* +1 )/ n,r )(0)}. We define 

n n (r,f) = Y f(g(^ i+1)/n ’ T) - 7 t ( * /ti ’ t) )). 
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From Proposition 2.5, we know .M>v( 0, r)flQ = 0. For each x € 0, r), 

set Z* = [nx\/n and = 1% + (1/n). Since there are exactly two dual paths 
7Tr X ’ r ) and 7r| x ’ r ^ starting from (x,r) with ^’^(O) > 7r[ x ’ T ^(0), from Proposi¬ 
tion 3.2(e) of Sun and Swart [29] it follows that : n € N} and ; 

n € N} converge to 7r|' t,T ^ and 7ri x,T \ respectively, in (II, dg) as n —> oo. 
Hence, 7 Z n (r, f ) —>• k(t, /) almost surely as n —>• oo. For each i G / n , there ex¬ 
ist Vi G (vr(V^)(o),7r((* +1 )/ n - T )(0)) nQ and G W such that G 

A4w(0, t). Hence, for uif = sup{|/(f)| :t> 0} we have 7 Z n (r, f ) < m/^w(0, r) 
for all n. As E[£w(0, r)] < oo, the family {TZ n (r,f) : n G N} is uniformly in- 
tegrable, and hence we have lim n _ K>0 E[7£ n (T, /)] =E[re(r,/)]. From the fact 

that ^( 7 r(( l + 1 )/ n > r )) _^( 7 j-(*/ n > r )) = IL(l/n + \/2IF) where IF denotes the stan¬ 
dard Brownian motion on [0, oo), we have 

lim E[7£ n (r, /)] 


= lim nE f(H(l/n + V2W))\l/n + min V2W{t) > 0 

n-*o o 1 fe[0,r] 


xpfl/n+ min -\/2IF(t) > 0^ 
V te[o,r] / 


= lim E 

n—>■ oo 


f(H(l/n + \/2IF))| min \/2IF(f) >-1/n 

te[0,r] 


n 


X (2$(l/v / 2rra) - 1) 

= E(f(V2W + ’ T ))/V^, 


where the last equality follows from Lemma 4.1, Slutsky’s theorem and con¬ 
tinuous mapping theorem. This completes the proof. □ 


Now, to complete the proof of Theorem 1.3 we need the following lemmas. 


Lemma 4.5. 
oo. 


For t > 0, we have T ) ra°’ 0 ^|l{L(o,o)>nr} =>■ V / 2IF +,T 


as n —>■ 


Proof. Using translation invariance of our model, we have 

- *(/(«. 

This holds for all / G C&[0,oo) which completes the proof. □ 


Lemma 4.6. For r > 0, we have: 

(a) sup{|T>i 0,0) (s) - T>i°’ 0) (s)| : s > 0}|l{ L ( Oj o)>nr} 40 as oo, 








28 


R. ROY, K. SAHA AND A. SARKAR 


(b) sup{|i4°’ 0) (s) -pD^ ,0 \s)\ : s> 0}|l{ L(0> o)> nT } 4-0 asn^-oo. 

Proof. For part (a), fix 0 < a < 1/2, T > 0 and we observe that 
P(sup{|Z)fc(0,0) — (0,0)| : k > 0} > n a ,L( 0,0) > nr) 

< P(max{|P)fc(0,0) — D/,( 0,0)| : 0 < k < n{r + T) + 1} > n a , 

L( 0,0) > nr) + P(L(0,0) > n(r + T)). 

Because of Theorem 1.2, it is enough to show that y / nP(nrax{|Hfc(0,0) — 
Dk( 0,0)| : 0 < k < n(r + T) + 1} > n a , L( 0,0) > nr) -A 0 as n —» oo. Here, we 
present the simple idea behind the proof; the details are available in Roy, 
Saha and Sarkar [27]. 

The distance d l k between 4(0,0) and the closest open vertex to the left of 
4(0,0) being n a or more has a probability (1 — p) n . Thus, the probability 
that the maximum such difference for 0< fc<n(r + T) + l is bigger that n" 
is of the order n(l — p) n . Similarly, for the distance d r k associated with the 
vertex rk (0,0). Since | Dk (0,0) — Dk (0,0)| <d l k + d k , as n —^ oo, we have that 
y / nP(max{|Z)fc(0,0) — £4(0,0)| : 0 < k < n(r + T) + 1} > n“,L(0,0) > nr) 
converges to 0. 

For part (b) of the lemma, we need T ) i°’°'|l{L(o,o)>nr} =>■ \/2 W +,T as 
n —>• oo which follows from part (a) and Lemma 4.5. Hence, r^(0,0) — 4(0,0) 
is of the order y/n. Also given 4(0,0) and rfc(0,0), the number of open ver¬ 
tices lying between these vertices has a binomial distribution with param¬ 
eters (rk(0, 0) — 4(0,0) — 1) and p. Since these open vertices together with 
4(0,0) and rfc(0,0) constitute (4(0,0), the proof follows from similar order 
comparisons as done in (a). □ 

Proof of Theorem 1.3. We remarked that VP 1 |[o,i] = fR +,1 |[o,i] = 
W + . The proof of Theorem 1.3 follows from Lemmas 4.5 and 4.6 and Slut¬ 
sky’s theorem with the choice of r = 1. □ 

4.2. Proof of Theorem l.f. For A > 0, let A := A 3 / 2 (\/27op) - 1 - We show 
that: 

Lemma 4.7. For t, A > 0, 

( OO 

L(0,0) > nr, ffCk(0, 0) > (An) 3 / 2 

k= 0 

l / r°° 

= — =f(V 2 W + ’ T (t)dt> A 

7ov 7rr V Jo 

= -- / F r+ (Xt- 3 / 2 )t~ 3/2 dt. 

270 V* Jr 0 
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PROOF. For / e <7[0, oo) let /(/) := f 0 °° H(f)(t) dt. Since P (W T € A) = 1 
where A is defined as in (21), I is almost surely continuous under the measure 
induced by W T on (7(0,00). The proof follows from Theorem 1.3(h) and the 
continuous mapping theorem. □ 

From the previous lemma, we derive the following. 


Corollary 4.7.1. For A > 0, we have 

i roo _ 

lim ##P(#<7(0,0) > (An) 3 / 2 ) =--= / F T+ (\t- 3 / 2 )t~ 3 / 2 dt. 

n ~>oo 27oV7T Jo 


Proof. For any r > 0, we have P(#(7(0,0) > (nA) 3 / 2 ) > P(L(0,0) > 
nr, #(7(0,0) > (nA) 3 / 2 ), and hence hminf^oo #nP(#(7(0,0) > (nA) 3 / 2 ) > 

WtJ JT D+ (\t- 3/2 )t~ 3/2 dt. 

We observe that 


y / nP(L(0,0) < nr, #(7(0,0) > (nA) 3/2 ) 


' L nT J 


<##P ^A.(0,0)>(nA) 3/2 


vfc=o 


< #nE 


LnrJ 


]T3 fc (o,o) 


fc =0 


(nA) 


-3/2 


= \Ar( [nrj + 1)E(5 0 (0,0))(nA)" 3 / 2 , 

where we have used the fact that {£#(0,0) = h k (r(0, 0))(1) — h k (l( 0,0))(1) : 
k > 0} is a martingale (see Proposition 2.3). From the earlier discussions, 
it also follows that E(Do(0,0)) < 2E((7) = 2(1 — p)p~ l where G is a ge¬ 
ometric random variable. Thus, limsup n _ >00 y / nP(L(0,0) < nr, #(7(0,0) > 
(nA) 3 / 2 ) = 0 as r —> 0, which completes the proof. □ 


Proof of Theorem 1.4. We first recall the result Lemma 6.1 of Resnick 
[25], page 174 which states that for nonnegative Radon measures //, /t n , n > 1, 
on [0,oo) d \ {0} we have fi n A fi if and only if (i n ({ 0, ici] x • • • x [0,a#) c —>• 
//([0,xi] x ■ • • x [0,Xd]) c for all x\,... ,Xd > 0 with (x\,... ,Xd) # 0. This re¬ 
sult implies that Lemma 4.7 together with Corollary 4.7.1 and Theorem 1.2 
prove (6). 

Fix r > 0, A > 0. For a < 2/3, 5 > 0 and for all large n, we have P(L(0,0) > 
nr, #(7(0,0) > (nA) 1 /") < P(L(0,0) > nr, #(7(0,0) > (n<f) 3 / 2 ). Fix any e > 0 
and choose b = 6(e) > 0 so that . )| ^ /7rr P(v / 2/ 0 ° o W +,T (t)dt > b) < e, where 
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6 = S 3//2 (joP ) 1 • From Lemma 4.7, we have 

limsupv / ralP(L(0,0) >nr, #C(0,0) > (nA) ly/ “) < e. 

n—>-oo 

On the other hand, from the properties of W + and VF T , it follows that 
P(/ 0 °° W +,T (t)dt > 0) = 1 for r > 0. Now for a > 2/3 and 5 > 0 we have 
P(L(0,0) > nr, #0(0,0) > (nA#/ a ) > P(L(0,0) > nr, #0(0,0) > ( nd) 3 / 2 ) for 
all large n. Again from Lemma 4.7, we have 

liminf \/nP(L(0,0) > nr, #0(0,0) > (nA) 1 /") 

n—>-oo 

> —\=F (V2 [ W + ’ T (t ) dt > 6) . 

V Jo ) 

Since 

limsup #nP(L(0,0) > nr, #0(0,0) > (nA) 1 /") 

< lim #^^(-^(0,0) > nr) = —^-=, 
n-*oo 7ov 7rr 

letting <5 —>• 0, we have limn^oo #nP(L(0,0) > nr, #0(0,0) > (nA) 1//a ) = 
for a > 2/3. This completes the proof of (8). 

The argument for (L(0,0), (ZD max (0,0)) 1 / 2 ) being similar is omitted. □ 
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